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Abstract 

We extend the method of differential renormalization to massive quantum field 
theories, treating in particular A(/? 4 -theory and QED. As in the massless case, the 
method proves to be simple and powerful, and we are able to find, in particular, compact 
explicit coordinate space expressions for the finite parts of two notably complicated 
diagrams, namely, the 2-loop 2-point function in Ay? 4 and the 1-loop vertex in QED. 
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Differential renormalization (DR)[1] is a coordinate space renormalization proce- 
dure which removes the divergences of bare amplitudes by writing these as derivatives 
of less singular functions and then prescribing the derivatives to be integrated by parts 
with all surface terms discarded. The singularities of amplitudes that did not allow 
for Fourier transformation into momentum space are eliminated and one ends up with 
finite, renormalized amplitudes that satisfy renormalization group equations. The stan- 
dard and simplest example of the procedure is for the 1-loop 4-point "bubble" diagram 
in massless A99 4 , where the following identity is used: 
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M is an integration constant which will become the subtraction scale of the renormal- 
ized amplitude. If we prescribe the laplacian to act only after integration by parts, then 
the r.h.s. above has a finite Fourier transform, and in fact defines the renormalized 
value of the (divergent) Fourier transform of the l.h.s.. Generally, one can say that 
differential renormalization provides a prescription to continue distributions defined al- 
most everywhere (i.e., bare amplitudes, undefined at a finite number of singular points) 
into bona fide distributions defined everywhere (renormalized amplitudes). 

The procedure exemplified above has been carried out at higher loops and in dif- 
ferent models [1,2,3]. One can furthermore check explicitly the form of the divergences 
being subtracted by examining the surface terms which contain them, thus verifying 
that the method indeed corresponds to a counterterm subtraction procedure [4]. 

In this Letter, we extend the method to allow for the inclusion of masses in an exact 
treatment. At the outset, it is clear that the presence of masses should not interfere with 
the method since the renormalization procedure is related to short-distance singularities 
whereas masses only change the long-distance behavior of correlators. In fact, one could 
simply take a pragmatic attitude and expand all massive propagators around zero mass 
and then proceed to use standard massless differential renormalization. However, the 
question remains whether massive field theories are amenable to a treatment exact in 
the mass parameter at each order in perturbation theory. We find that our expectations 
are entirely fulfilled. In analogy to the massless case, renormalization is accomplished 
through the use of (massive) differential renormalization identities and an integration 
by parts prescription. We treat QED at one loop and A<p 4 -theory at two loops and, 
most notably, we are able to give in closed form the full renormalized expressions for 



2 



the 2-loop "setting sun" diagram of A92 4 and the 1-loop vertex of QED, two results of 
considerable difficulty of calculation in standard momentum space treatments. In the 
end, we shall discuss the general features of our modus operandi. 

First of all, let us recall the euclidean propagator for a massive scalar particle in 
four dimensions, 

A ( .,m) = — — (2) 

where K\ is a modified Bessel function. Up to permutations of external legs the only 
1-loop diagram contributing to the 4-point function in \ip A theory is: 

( 4) _ A 2 / 1 mK^mxW 

At short distances, this has the same log divergence as the corresponding massless 
diagram; in the spirit of differential renormalization we look for an expression that 
corresponds to the massive generalization of Eq.(l). We find: 

m 2 K 2 (mx) 1. A n.mKo(mx)Ki(mx) 2n M 2 r( - 4 \ . . . A , 

^ L = -(a-Am 2 ) - '— ^ +TT 2 ln^ T S {4 '(x) (4) 
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where M = 2M/7 and 7 = 1.781072... is the Euler constant. The contact term has been 
added in order to give a well-defined massless limit to the r.h.s. above, coinciding in fact 
with Eq.(l). We shall further analyze this later on, and in fact we will see that Eq.(4) 
can be seen as a prototype in our treatment, since it presents all the basic guiding 
elements in working out massive DR identities in general. The renormalized 1-loop 
4-point function, T^\x, M), is gotten by simply substituting the above in Eq.(3). 

The one-loop beta function can be obtained from the renormalization group equa- 
tion 

{ M jM + ^T\ +lmm2 ^T 2 - Al ) [-^ (4) (*) + rg>(s,M)] = 0, (5) 

where /3(A) is the /3-function, 7™. (A) is the anomalous mass dimension, and 7(A) the 
anomalous dimension of if. Dropping terms of higher order in A (viz., j m and 7), the 
result for the /3-function is: 

« A >= 3 (ii) 2 < 6 > 

with the three permutations of the s, t and u channels having been added. Our result 
for the bubble diagram can be easily Fourier transformed. The result in momentum 
space is: 



A 2 
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(7) 
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which agrees with the standard result in textbooks [5]. Note the appearance of the 
2-particle threshold 4to 2 in the operator which was extracted in Eq.(4). This will be 
a recurring feature in what follows. We have also worked out this amplitude for two 
different masses running along the two lines in the loop; we do not present it here, 
but rather just mention that the result generalizes the one above and, as expected, the 
threshold (m\ + 771-2) 2 appears instead of 4m 2 . 

A far less trivial example showing the power of differential renormalization is given 
by the computation of the 2-loop correction to the 2-point function (often called the 
"setting sun" ) . The complexity of the finite parts of this diagram is such that they are 
not presented in standard reviews of \tp 4 renormalization. The bare expression for the 
diagram is: 

(2) , . _ A 2 / 1 mKximx] 



This bare amplitude is quadratically divergent and thus requires the extraction of two 
laplacians in order to have a good Fourier transform. Little effort is needed to verify 
the identity that leads to the following renormalized value of the above diagram: 

A 2 



rg } (x,M) = m ^2y ( a -9m 2 )( a -m 2 )(rn 2 Ko(mx)Kf(rnx)+m 2 K$(mx)^ 

(9) 
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This result is by itself a remarkable application of differential renormalization. It corre- 
sponds to a closed and compact expression for the renormalized 2-loop 2-point function. 
Its Fourier transform is complicated and we will not work it out here. Again, we note 
the appearance of the 3-particle production threshold as a healthy sign of the proce- 
dure. Further, the second operator takes the form (□ — to 2 ), which vanishes on mass 
shell. The coefficient of the DS A (x) term is fixed as in the previous example by the re- 
quirement of a smooth massless limit and this, in turn, fixes 7(A) in the renormalization 
group equation to its standard value 

7(A) = ^f^-V. (10) 
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On the other hand, the coefficient a of m 2 5^{x) is not fixed by a strict massless limit, 
which is to be expected, since it corresponds to a mass subtraction, and fixes the 
(scheme-dependent) "f m function in the renormalization group equation. In particular, 
a can be chosen to be —1 so that the whole 2-loop contribution vanishes on mass shell. 



4 



Though we do not intend to present an exhaustive computation of the perturbative 
expansion of massive Xip 4 , let us comment on a few more diagrams. There are only 
two more diagrams contributing at two loops to the 4-point function. The first one, 
usually called the "ice cream cone" , though extremely difficult in momentum space, is 
easily calculated along the lines we have sketched above. The second one, consisting of 
two 4-point 1-loop bubbles attached together (and thus, the "double bubble") is not as 
straightforward due to the presence of a convolution. 

We now turn to the computation of 1-loop diagrams in QED 0. The (massive) 
fermion and photon propagators are, respectively: 

„, . 1 .mK-\(mx) , . 

S( x ,m = — 2 11 

and 

this latter one being given in Feynman gauge. The ipJ^ij) vertex has the value ie^^. The 
1-loop fermion self-energy then reads: 

/ % / ie\ 2 I . mKi(mx) . . 

Z bare (x) = ^—j - m) 1 j— L ^ (13) 

The mass piece above has a logarithmic short-distance divergence, while the derivative 
term is linearly divergent at short distances. These divergences are eliminated by the 
following massive DR identity: 

~ 3 ^ = -(D-m 2 ) °\ J +7T 2 ln—5^(x), (14) 

where again the contact term is determined such as to give a well-defined massless 
limit to the r.h.s. (and again equal to the corresponding massless identity). The final, 
renormalized expression for the fermion self-energy will then be: 



K (mx) m 2 tfK (mx) 
[D — m ) \{a + 4m) 5 h — 5 



e 2 
32vr 2 l 

+2tt 2 \ii -%{^ + bm)5^{x) 



x' 2 



(15) 



The non-local piece of this amplitude vanishes identically on mass-shell due to the 

operator (□ — m 2 ), as expected. In analogy to the previous case, the $5 4 (x) contact 

"'"For a complete set of momentum space techniques to compute one-loop diagrams in gauge theories, 
see [6]. 
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term is fixed by the massless limit (and leads to the standard value for the anomalous 
dimension of the fermionic field), whereas the m5 4 (x) term is not. We leave the coeffi- 
cient b undetermined, and again it is clear that a particular choice, viz. b = — 1, makes 
the whole amplitude vanish on mass-shell. The Fourier-transformed, momentum space 
amplitude can be obtained fairly straightforwardly, and we do not present it here. 
The bare 1-loop vacuum polarization is: 



.mK\(mx) , „ . mKi(mx) 
m) -Ju(-^ ~ m) - 



(16) 



x x 
By standard manipulations with Bessel functions and the use of a massive DR identity, 
Eq. (4), we find the following renormalized vacuum polarization at one loop: 



U Rfll/ (x,M) 



24vr 4 



+—(K 2 (mx) - K 2 (mx)) ) + 2tt 2 In ^j8^(x) 
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(17) 
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As expected, the result is automatically transverse and shows the presence of the 2- 
particle production threshold in the operator coming in front of the non-local piece of 
the amplitude. 

We finally turn to the 1-loop vertex. Its bare value is: 

1 



fe / ie \ 3 mKi(mx) mK x (my) 
( x > y > = {~^2) "?pW- m > 7/i (-$-m) 7/ , 



x y ■ (x - y) 2 ' 

In the above, the only divergent piece is the one containing the two derivatives: 

1 



Vf v (x,y) 
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We renormalize this by an identical procedure used in the massless case [1], i.e., by 
integrating the two derivatives by parts onto l/(x — y) 2 , and separating that into a 
(finite) traceless piece and a (divergent) trace piece. The divergence, thus isolated in 
the trace, is then renormalized with massive DR identity Eq.(4). The final result reads: 



V?(x,y,M) 



( ie \ 3 
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(x — y) 2 \dy^ dx^ 2 
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This concludes our results. We have extended the method of differential renormal- 
ization to massive \ip 4 and QED, and the results we present here, though not intended 
to be exhaustive, should indicate both the feasibility of, and the guidelines for, the 
treatment of massive fields at higher loops and in different models. These guidelines 
basically are: a) the presence of an equal number of K-functions on both sides of 
massive DR identities; b) a direct connection to massless differential renormalization, 
through the substitution of lnx 2 M 2 in this latter case for Ko(mx); c) the appearance 
of the appropriate particle production thresholds in the differential operators present in 
massive DR identities, and d) the appearance of contact terms, due to the requirement 
of existence of the massless limit for renormalized distributions. These features have 
been of central importance in helping us find all the renormalized amplitudes presented 
here. 

Two final comments are in order regarding contact terms in massive DR identi- 
ties. Firstly, the presence of these terms may seem puzzling, since they vanish at 
points where massive DR identities are valid mathematical equations, i.e., for i / 0. 
However, insofar as these identities are prescriptions for extending (i.e., renormalizing) 
distributions, they should be understood as identities everywhere, including at contact, 
and specific contact terms will then determine precisely and uniquely the value of the 
extended (i.e., renormalized) distributions at those points. Finally, it is also worth 
noting that contact terms encode the entire renormalization group freedom present in 
the renormalization procedure. Our choice has been to fix contact terms in order to 
agree with the massless case, which corresponds to a wave function renormalization 
prescription, and to leave undetermined the mass subtraction, though an on-mass-shell 
scheme has been indicated. 
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